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ABSTRACT 


Real time simulation of the time-dependent behavior of low power 
thermal reactors is amenable to analog computer techniques, This paper 
presents the basic physical concepts relative to thermal reactor operation, 
a derivation of the governing transient equations noting simplifying 
assumptions, and the development of the analogous electrical circuit to 
ultimately provide solutions to these equations, In addition, many —Ó 
problems involved in the simulation of & specific reactor are discussed 
and analyzed in detail, As a result, the successful design and con- 
struction of a simulator whose time-dependent behavior closely approxi- 
mates that of an actual reactor is demonstrated. The reactor simulated 
is the Aerojet -General Nucleonics low power thermal reactor (AGN-201) 
currently installed at the United States Naval Postgraduate School, 
Monterey, California, Flux-time recorder plots for various sequences of 
operation of reactor and simulator are included to indicate the degree of 


simulation attained, 
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A LOW FOWER REACTOR SIMULATOR 


1, Introduction, 

Detailed analysis of the transient behavior of nuclear reactors 
has been greatly facilitated by the use of electrical or electronic 
simulators, A reactor simulator is essentially an analog computer with 
circuitry specifically designed to solve, on a true time base, the 
rather complex differential equations governing time-dependent reactor 
operation, 

As will be evident from the development of the governing equations, 
complete solution by direct analysis would be extremely laborious, if 
not impossible, for other than a few limited cases where modifying 
assumptions must be made before a useful solution can be obtained, 
Moreover, direct mathematical analysis could not be expected to instill 
in one any real feeling for reactor kinetics under typical operating 
conditions, Because circuit values can be measured quite precisely, 
accurate and continuous solutions to the transient equations for almost 
any imposed operating conditions may be readily obtained from the sim- 
ulator, In addition, by employing a simple simulator control system, 

a strong feeling for reactor responses to any of a multitude of typical 
Operating situations may be achieved, This particular feature of the 
simulator might permit reactor operating personnel or students to receive 
the equivalent of a great deal of reactor control experience and 
operational training when access to an actual reactor is not feasible, 

Because analog computers are well-tried and accepted engineering 
tools for use by those involved in the study and design of dynamic systems 
of all types, reactor simulators find application in preliminary reactor 
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design and especially іп the evaluation and perfection of associated 
systems which might be ultimately employed to control an actual reactor, 
In other words, completed control components and systems may be connected 
to the simulator by means of electronic transducers thereby permitting 
operational checks to be made on their specific design functions prior 

to actual installation, 

Depending on the complexity of the system to be studied and the 
degree of simulation desired, simulators may range іп size from small 
units such as the one discussed herein to extremely large, complex, and 
expensive systems involving the use of great numbers of amplifiers, function 
generators, servo and electronic multipliers, electro-mechanical devices, 
and countless precision components, In the case at hand of a particular 
low power thermal reactor, utilization of the basic simulator circuit 
together with a relatively simple control network simulates actual reactor 
operation to a surprisingly good degree of accuracy, 

The discussion of basic concepts, the derivation of governing 
equations and the build-up of the electrical analog represent a compilation 
of treatments given on this subject by the first seven authors listed 
in the bibliography. The simulator itself is essentially a Pagels 115 
simulator, modified to duplicate more closely the responses of the 
individual reactor under study. 

It must be remembered that the degree of simulation attained can be 
no better than the degree to which the dífferential equations utilized 
define the system under study, The transient equations employed in 


thermal reactor simulation contain many simplifying assumptions and 


Numbers in brackets refer to Bibliography in Section 8, 





approximations; and as a consequence, precise duplication of actual 
reactor behavior may not be expected, Despite these assumptions, however, 
it will be seen that simulator solutions to the equations employed do 
closely approximate observed reactor behavior, 

It is the purpose of this thesis, therefore, to present a rather 
complete study in the design of a reactor simulator with the ultimate 
goal being to construct such a device and then compare time dependent 
behavior of the actual with that of the analogous system, 

2. Basic Thermal Reactor Concepts, 

In order that the discussion of the simulation of a low power 
thermal reactor be complete, it is important that basic concepts be 
briefly reviewed prior to the derivation of the governing differential 
equations, 

When the nucleus of a fuel atom fissions, it splits into two primary 
fission fragments of more or less equal mass accompanied by the prompt 
emission of from two to three neutrons of high energy. On the average, 
about 2,5 neutrons are liberated per thermal neutron absorbed in físsion 
of a 0232 nucleus, Тһе majority of these have епегріев of about one or 
two Mev, Most of the fission fragments thus formed are radioactive and 
some are themselves high energy neutron emitters, The decay constants 
for six delayed neutron precursors have been measured although the precise 


origin of only three of the six is currently known, namely, BRÊ” 


„ BR”, and 
m Ihe average energy of the delayed neutrons is about ,5 Mev, The 
properties of the delayed neutrons play a vital role in time-dependent 
reactor behavior and permit safe control of nuclear reactor fission 


processes which might otherwise be extremely difficult, 


In addition to having small thermal neutron absorption cross-sections, 





good moderator materials are those made up of elements of low mass number, 
where the energy transfer from neutron to moderator nucleus on collision 
may approach 100% as in the case of a hydrogenous moderating medium, 

This results in extremely rapid slowing down of the high energy neutrons 

to so-called thermal energies (about ,025 ev), At this point, the neutrons 
ane essentially in thermal equilibrium with the moderator atoms (or mole- 
cules) and hence are referred to as thermal neutrons, Despite the 
extremely short slowing down time, some of the neutrons will either 

escape from the reactor or be absorbed by the core material before 
thermalization so as to be unavailable for subsequent fissioning or neutron 
regeneration, 

Once the fast neutrons have been reduced to thermal energies, they 
may leak or diffuse out of the reactor or be absorbed by moderator, fuel, 
or poison nucleii, A large fraction, about 854$, of the thermal neutrons 
absorbed by um 0232 fuel will cause fission; and a new neutron lifetime 
cycle begins again, The next job, therefore, will be to express in equation 
form, the lives of neutrons as described above, 

3, Derivation of the Transient Equations, 

Based on the conservation or balance of neutrons principle, applicable 

to any reactor, a general expression for the rate of change of neutron 


population density (n) with time can be written as follows; 


Production-Absorption-Leakage = 2. (1) 


The rate of production of neutrons may œ obtained from empirical con- 
siderations, Likewise, the rate of absorption depends on various absorp- 
tion cross-sections which must be determined experimentally, For determin- 
ing leakage, resort is made to classical mechanics and the familiar laws 


of diffusion theory which, although somewhat more complicated, can be 


4 





reduced to simplified expressions when certain prudent assumptions are 
introduced, Let us now form the basis on which the various simplifying 
assumptions can be made, 

The thermal reactor to be considered here is of the homogeneous type 
with reflector, The fuel is a homogeneous mixture of 053” in the form of 20% 
enriched U0, with polyethylene, (сњ), as the moderating material, A thin 
high density graphite reflector is employed, However, in the analysis 
which follows, it will be assumed that the reactor is homogeneous but un- 
reflected, The degree to which this assumption is valid can best be 
ТИН by a comparison between final simulator solution and actual 
time-dependent reactor test runs, As will be seen later, the assumption 
was quite reasonable in this case, 

Before continuing with additional assumptions, it will be helpful to 
rewrite equation (1) as the familiar diffusion equation for thermal 


neutrons in non-equilibrium, basic to reactor theory, namely 
| on 
am 9 መ = 2 = 
5 25,9 (-рУ 29) DE (2) 


where S is the production term in neutrons per cm? per second, E a is the 
macroscopic absorption cross-section for monoenergetic (or thermal) neutrons 

in “ай $ is the thermal neutron flux in neutrons per cm? per sec, and 

(-DV 2p) is the neutron leakage per cm? per sec, By the acceptance of 

equation (2) as the equation representing (1), it has been assumed that we are 
dealing only with monoenergetic (or thermal) neutrons and then only at dis- 
tances greater than two or three mean free paths from strong absorbers, 

sources, or boundaries, This also assumes that the moderating or diffusion 
medium is a poor absorber relative to its ability to scatter and, hence, slow 
down the high energy neutrons, The partial derivative of n is utilized because 


n is a function of both position and time, However, it will be presumed 





that the reactor is at or near critical and the shape of the flux 
distribution remains constant even though it changes in amplitude with 
time, Under the latter provision, (-DV?$) could be represented by some 
constant multiplied by the flux $, Making use of this, therefore, let 
(-руаф) = 2 ¢, where the symbol Z will be used throughout this 
development 5 indicate the macroscopic cross-section for or the probabil- 
lty of occurrence of the particular process indicated by the subscript, 
i.e., in this case - 1 to mean probability of loss by leakage or leakage 


coefficient, This essentially eliminates the space variable and permits 


dn 
dt 


under the above considerations is merely the product of n, the density 


the use of the total derivative of n, or in lieu of = The flux ф 
of thermal neutrons in neutrons per cubic centimeter, and v, the mean 


thermal neutron "speed." Or, expressed in equation form, 
® = nv neutrons /cm?/sec 


Thermal neutron flux $, therefore, is a scalar quantity, We are now in 
a position to examine each of the terms of equation (2) in an attempt to 
develop a more workable and meaningful equation, 

The source term (S) must be broken into several parts, i.e., Sn S4» 
and Sg, referring respectively to prompt, delayed, and those uniformly 
distributed as & consequence of an extraneous outside source (e g., a 
Ra-Be source), Expressed in equation form: 

5 5 5, «5, %% | (3) 

The rate of fission can be expressed as ይቃ. If V represents the total 
fission yield of fast neutrons from both delayed and prompt sources and 


p represents the total fraction of these neutrons which are delayed, then 
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the prompt neutron generation rate 1s given by 


2,9 Y (1-8) 


However, not all of these neutrons reach thermal energies, As previously 
mentioned some leak out or escape and others may be captured by fuel or 
other nucleii during the moderating or slowing down process, Probabilities 
have been evaluated for each of these processes, The former is given by 

the fast non-leakage probability, qur (where B? is the so-called 
buckling factor related to size and T is the Fermi Age, a measure of the 
distance neutrons travel in slowing from fission to thermal energies), and 
the latter by p, the resonance escape probability, For simplicity, these 
factors will be combined into a single thermalization probability symbolized 


by Pth’ Therefore, the prompt thermal source becomes 


5 = фа, (1-8) Vea (4) 


-p2 
where Pin = ре 3 E. Now the remainder of the neutrons are the so- 


called delayed neutrons associated with the particular fission fragments 
mentioned in section 2, The delayed neutron generation rate, therefore, 


ls given by 


Zo gv 


Now, as previously stated, six delayed neutron precursors have been 
determined so that the delayed neutron generation rate of each of the six 


may be expressed by 


е: Хор 


where В; = fraction of total number of neutrons contained in the ብ 
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delayed group, It should be noted here that в = 2 At any time, 
Tz 


there will be a concentration C; of the шча group of delayed neutron 


emitters present in a unit volume of the core, Because these delayed 


neutrons are emitted as a result of the normal radioactive decay of the 


precursors, each emitter obeys the familiar radioactive decay law, In 
th 


other words, the rate of change of the concentration of the i group will 
be equal to the production rate minus the decay rate, or 

dc; 

ae cow WWE 
where N is the decay constant for the а emitter, 
Therefore, 

ас 

۳ 
dt ` Бэк - 5,6, (5) 


Delayed neutrons are therefore supplied to the reactor at the same rate at 


which decay of the precursors takes place, This can be expressed as a 
6 
summation 3 A, Oy. Because the average energies of the delayed 


1-1 
neutrons are somewhat less than those of the prompt, there will be less 
tendency,for the delayed neutrons to leak out during the slowing down 
process; and hence they will have larger values of non-leakage probabilities, 
The difference in the: values of non-leakage probability for the two types 
may be taken into account by utilizing effective values of B, in the 
final equations which are slightly larger than the actual values, This, 
then, permits multiplication of the delayed neutron term by the same value 
of non-leakage probability as that used for the prompt neutron term, 


Combining this factor again with the resonance escape probability p, 





the rate at which delayed neutrons are supplied can now be expressed as 


a summation of terms, namely 


6 
it hc uu №; (6) 
121 
where 4,0, із defined by equation (5). Finally, бо comes from an 
essentially constant extraneous source of neutrons which makes a contribution 
to production rate in the form of a uniformly distributed source, This 
term will be discussed in more detail later, 

The absorption term in equation (2) needs no additional explanation 
except to say that 24 represents the macroscopic absorption cross- | 
section for all processes, that is absorption of thermal neutrons by 
structure, moderator, fuel, or poison nucleii, 

The leakage term was discussed earlier and can be taken to equal 

Ра 19, where » 1 18 dependent on the moderator characteristic and 


reactor size, 


The component terms of equation (2) can now be assembled as follows; 


dn _ 
dt E (5, + За + 56 ) - 2. ф =- 2. Ф 


d = 
ችን ә 2,01 - B) V Pan б Pth > XC, * ы - ф 24 P 20 
i=] 


or 48 bz, ]1 - ۷ ር. DUM - ቁርጆ, ፥2)) + & (7) 


ዘ 


In order to reduce equations (5) and (7) to more usable forms, several 


terms must be defined, First, the term reactivity ሥ may be defined as 


К 
eff 
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where K er is the effective multiplication factor for a given finite sized 
reactor, By effective multiplication factor is meant the population ratio 
of successive neutron generations, If 1” is next defined as the mean lifetime 
of thermal neutrons in the reactor and 1/1, as the rate at which neutrons, 
prompt and delayed, are thermalized in the reactor, then К ер could also 


`+ 
be expressed as the ratio of 1 p і 6 / , or 


K = 1A, (9) 


This may be interpreted as follows: Kore 1s the ratio of the time spent 

by the thermal neutrons diffusing in the reactor to the time spent in being 
thermalized, In other words, if it takes longer to thermalize the second 
generation of thermal neutrons than the life expectancy of the first 
generation, the chain can be expected to die out, In this case, К Pr would 


be less than unity, The reverse would be true if K were greater than 


ff 
unity, When K pp 1s precisely unity (or in the case of the presence of 
an extraneous source, slightly less than unity), just as many thermal 
neutrons are born as die and the chain 1s sustained, These three situations 


represent subcritical, supercritical, and critical behavior of the reactor, 


er can be expressed as the reciprocal of the thermalization rate [6] , "UT 


№ = ea (10) 
У.У ву, 

This term is also given the name of modified lifetime by Murray Ex 2 

Now thermal neutron lifetime іп ап infinite medium would be given by the 

average distance traveled by the neutron before absorption, called the 

absorption mean free path, № ( ው divided by the mean neutron 


velocity v. However, in a finite reactor, this must be reduced by the 
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thermal leakage factor 


Thermal Absorption 
Thermal Absorption + Thermal Leakage 


ዙ፦ . .( _Û 1 р Za ж 1 
BE 7-4 эссе ፆ* 


Substituting (9), (10), and (11) into (8), therefore, 


— 1 ۷۱2 2 й 1) е 
— 1 ጩጨ ۱۳ я መ 1 em ". 
/ QUAE ۲ م2‎ ۷ Pin ся 


Now, utilizing 1", 1, and JP as defined above, equations (5) and (7) may 


be simplified as follows; 








6 
dn _ "۳ 
EU = lez, $1 - gv + я UN ? (ይ. +2) * dà 
dn 1 3 д 
ап - n (1 - 8) +р C ar 05 
dt ME th 6 14 5 
п 1 
= (ТЄ в - — | +p E С + 5 
ds 6 
9 = P T Ean ее "5 (13) 
ac 
Likewise, from (5), ون و‎ Б (=, Әр, = i ) 
P, ас, 
ог Мб, Se” Im ማሺ 
lo P 7 
th 


Equations (13) and (14) are the desired transient equations, 


ጊጊ 





Finally, a single governing equation can be written by a simple 


substitution of (14) into (13) as follows; 


ас 
і 
d 2 5 ዞ4 
E: £— з, г ux 3867... 
M lo p 
ms th 
в. ЕТ Ме oo 
las 7. መክ 5 "ክፍ £e ар * “е 
1o ዕ i= 
چم‎ * = 
= D 7 Pth шо" Зб 
15 = at 
ас 
* dn п x 6 m x 


It will be further assumed that the so-called one-group theory 5 
sufficiently accurate for accomplishing the goal at hand, In one group 
theory it is supposed that all production, diffusion, and absorption of 
neutrons occur at a single energy, namely, the thermal energy, Although 
this is a considerable approximation, the final simulator circuit is 
thereby simplified and results are still quite satisfactory, Since this 
theory implies that all fission neutrons have thermal energies, the 
resonance escape probability, p, and the non-leakage probability during 
the slowing down process, «ВТ » do not affect the multiplication factor 
and, hence, become unity in equations (13) thru (15), The two terms are 
carried through in the development of the electrical analog, but no 
attempt is made to employ them as variables in the final circuit, 

4. (Тһе Electrical Analog. 
Equation (15) is now in a convenient form for a computer solution, 


The network given in Figure 1 will be shown to describe this differential 
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equation, Referring to Figure 1, and using Kirchoff's second law, it 


is evident that 


6 


I, З & - б ود‎ + I 


Note the similarity in form to (15), 


Is ED 
The Rate Term `© dt 


and letting E in volts be proportional to n in neutrons per cm? per 


sec by the constant Y و‎ 1,0. Е = ыў n, then 
с + 
І, = n: „ ፳፻ (16) 
lo 


Тһе Reactivity Term (on) 


Let A be the proportional distance of the potentiometer tap (at P) 
from the center of the potentiometer, Referring to the more detailed 


development given in Appendix I, 


AB 
8 R ፡ + R/P + PAR (1 - a2) | 


The range of values of A are zero at the center, +1 at the +E end, 
and -lat the -E end, Now, it is apparent that I, will have its 
maximum positive value at A = + 1, Therefore, the range of reactivity 


variation of the simulator may be set and the proper scaling factor 


determined by setting 
14 





max 8 пах 
0 ۵ a 
or ү ፊስ pi Га жы (17) 


* 4R(R + P) 


Since I, max occurs at A = 1, then 


E 


I, max 5 В(1 + 


and therefore, 
۳ E 
МИ кент — en 
/ max 


Again representing E by J n, this reduces to 


ብ | መሥ" “= 


More careful consideration of (17) above reveals е is a non-linear 
function of A which is undesirable for our purposes, This is due to 

the potentiometer loading having its greatest effect around the mid-position 
of the potentiometer, Linearity can be very closely arproximated, how- 
ever, by making R sufficiently large in relation to P that the 


denominator of (17) becomes essentially unity, leaving 


ሆው. у= — 
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If, for example, R=5P, the linearity remains below one percent for any 


value of A from (-)1 to (+)1, 


6 
ЗР i 


The Delayed Precursor Term 


6 


` The ж term of (15) is simulated by the six Rc, delay networks, 
1- 





The basic voltage equation describing the voltage drop across any of the 


n type networks is 


E - Ww + R а, 
5 с, 47 dt 


where 9, is the capacitor charge in coulombs and Cy is the capacitance 
of the на group. The electrical charge is analogous to the "storing" 
of delayed neutrons in the fission fragments before they are emitted, The 
corresponding delayed neutron equation represented by these networks is 
(14), which upon rearranging becomes 

10 Pin X6, lo Py, 40 


1 
n = سس‎ MÀ— MÀ лы + — — 
By Ва аў 


Again letting E =3n and matching terms with the basic equation above, 


4 35 
o Pin 1 
dQ J 
i ж ас 
Now, I, = ጭ፦ = (—]T—) lp, i. 
ас, 
and by substitution for iv. from (14), 
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J 


BR, 


But, ть 
6 
۱ 4 x > Е Б 2 x 
and Ны 1% 545. ащ. > ) ER | В; n = EA 455. М ча (19) 


Therefore, B, R, u = BR - Б, ..... 


5 + 
1 = 1 ) в, є- CR) do Pin МО 


The Source Term (б...) 
i ከ ቃሪ Ро, then 


AKT “а ኝ 
Ig rn um 
where Ag is the fraction off ground potential and V is a constant 
voltage source, Now, letting Ag = = then 
% (пах) 
ጁ 
Io - я lo Sg (20) 
Қо Зо (пах) 0 


Finally, summing the currents as given by equations (16), (18), (19) 


and (20), 
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тьма INC A Dp al Pata Apte 


her 
و(‎ 


— ا‎ 4, (21) 
Ro9o (max) lo 


In order for the current output of each term to be of the correct 


magnitude, the scaling factors which are enclosed by brackets in 





b equation (21) must be identically equal to each other; that is, 
3 
3 
“А 
& ml. З - об 1-3 й ( (22) 
e EA ГЕ” 52 5 
S| 1, RU + Р) пах ВЕ, ROS (max) lo 
жщ 


If this can be done, then these coefficients drop out of equation (21) 
leaving equation (15) which is to be duplicated electrically, Values to 
insert into (22) will be determined by constants in the actual system to 
be simulated as well as arbitrary values selected in accordance with standard 
computer operating techniques, 
5% Simulation of a Specific Low Power Thermal Reactor —AGN-201 

General thermal reactor simulation, therefore, can be readily duplicated 
by the circuit given in Figure 1, This circuit when set up on a conven- 
tional analog computer, using completely arbitrary circuit parameters, will 
yleld typical time-flux curves, The problem at hand, however, is to 
simulate a specificreactor wherein lie many additional considerations which 
first must be resolved before a final circuit can be drawn up, These 
include detailed investigations to, 

(1) determine actual range of reactivity variation, compute rates 
of reactivity change for control rod movement, and provide suitable means 


for duvlicating control rod action; 
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(2) determine a means of introducing a maximum reactivity control 


and temperature coefficient of reactivity;. 


(3) solve electrical problems resulting from real time simulation over 


wide operating ranges of neutron flux from source level to full power; 


(4) and to determine an appropriate value of neutron source voltage 


corresponding to the extraneous source production rate, 


In the analysis of these problems, attention shall be restricted to 


the simulation of the low power thermalreactor currently installed at the 


United States Naval Postgraduate School, 


referred to as the AGN-201, 


This reactor will hereinafter be 


Before solving the problems, it is necessary to mention a few essential 


features of the AGN-201 which will form the basis on which ultimate 


simulator construction may proceed, 


reactor are listed in Table 1 below, 


The general characteristics of the 201 


TABLE 1, GENERAL CHARACTERISTICS AGN-201 


Power (licensed) 
Maximum Thermal Flux e 100 mw 
Critical Mass 
Core Size (cylinder) 
Moderator (polyethylene) 
Reflector (graphite) 
Shield (lead) 

(water) 
Over-all diameter 
Over-all height 


Temperature coefficient of reactivity 


(based on water shield temp, data) 
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100 miliwatts 

4.5 x 10°em™? sec”! 

665 gm 023? 

10" dia, x 94" high 

11 kgm 

20 cm thick (1,70 gm/em?) 
10 em thick 

55 em thick (1000 gallons) 
6} ft 

9 ft 

-2.75 x 107^ ዶር 





Figure 2 shows the AGN-201 core tank containing the uranium-polyethylene 
core and portions of the graphite reflector. Figure 3 is a cross-sectional 
view of the entire reactor unit showing the core, reflector, lead and 
water shields, glory hole, and control rod assemblies, The critical mass 
consists of about 645 grams of y23? in the form of 20%- enriched U0,, It 
is contained in the form of a powder dispersed in a solid polyethylene 
moderator. The core is made in nine sections of various thicknesses 
(Figure 2). Each of the four bottom discs has three 5-cm holes and one 
2,5-cm hole, The 201 has two safety rods and two control rods, The 

three larger holes are for two safety rods and one coarse control rod while 
the smaller hole is for a fine control rod, The two safety rods and the 
coarse control rod are identical and each contain about l4 grams of 0232, 
The fine rod containsabout 3,5 grams of 0232. All rods operate so that 
reactivity is increased as the rod is inserted, The amount of reactivity 
each rod controls increases as the amount of active material contained 
increases, With the above amounts of 023? in the rods, the safety and 
coarse control rods control about 1,1% reactivity while the fine control 
rod controls about 0,3% reactivity, The rods are inserted by a reversible 
d-c motor drive mechanism requiring about one minute to travel the full 
25-сп stroke, Referring to Table ول‎ the thermal flux at licensed full 
power is given as 4.5 x 10° neutrons /cm?/sec, This is well below that 
at which Xenon poisoning becomes a consideration, namely, at a flux in 


12 neutrons /cm?/sec, Finally, because of the low 


excess of about 10 
power developed, heating of the core is so slight that it is essentially 
ambient temperature which determines core temperature and, hence, the 


maximum reactivity available, The negative temperature coefficient listed 


in Table 1, therefore, is based on the change in temperature of the 1000 
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gallon water shield surrounding the core tank ўзо! . With these features 
in mind, let us proceed with the investigations into the special problems 
mentioned, 

First, a control rod calibration must be made, the result of AS / / 
shown graphically in Figure 4, The term relative reactivity as used here 
implies an arbitrary zero reactivity reference when the rods are fully 
retracted, As will be seen later, the major portion of the reactivity 
contribution or "worth" of the rods is employed in reducing the amount 
by which the reactor is subcritical, Only a small fraction of this 
relative reactivity is actually positive reactivity as defined by (8). 

This curve represents & superposition of data taken on both fine and 

coarse control rods which actually exíst on the reactor, thus providing 

a single curve to represent the movement of both rods simultaneously. 

In view of the fact that the action of the two control rods on the reactor 

can be duplicated more simply and more economically on the simulator by 

a single controlling device, simulation of operation of each rod individ- 

ually does not seem warranted, Now, under normal conditions, the reactor 

will be operating in the straight line portion of the curve when near- | 
critical conditions exist, As indicated in Figure 4, the slope of thís 

straight line 1s 237 p, (ox 1079) per cm of rod travel, From weight 

data on the fuel discs which make up the two safety rods and the fine 4 y) 
and coarse control rods, it is found that the relative reactivity con- У 
tribution of each safety rod is essentially equal to that of the coarse S 


N 


N 

control rod or about 11,275 yo each, Therefore, the total contribution 4 У 
Ў N 

In terms of relative reactivity of all rods is found to be 36,725 ME : м 


The reactivity potentiometer can now be designed to match the control and . 


safety rod curves, Figure 5 shows a plot of percent of total resistance 
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of the reactivity potentiometer versus percent of total rotation (about 
350°) where all rods are being considered. Notice that the curve is 

shown up.to the point at which the safety rods are completely inserted 

as a dashed line. Because the safety rods are not utilized for control 
purposes, their individual calibration curves are of no concern as long 

as their ultimate contribution to total relative reactivity is known, 
namely, 22,550 МО ог about 61.4% of the total available for addition 

by rod insertion. From this point on, however, it would be desirable to 
follow the remainder of the curve which is essentially Figure 4 super- 
imposed but plotted using non-dimensional relationships. To duplicate the 
entire curve precisely would require a non-linear potentiometer which could 
be designed but at considerable expense. Instead, in view of the fact that 
the straight line portion of the combined control and calibration curve 

is of greatest interest, a center-tap potentiometer permits precise duplica- 
tion in this region. At 50% rotation, extrapolation of the straight line | 
portion back to this point reveals that 25% total resistance should result 
at this point. For the 20,000 ohm potentiometer selected, this requires a 
5000 ohm resistor in parallel with the first half of the potentiometer. The 
arrangement is shown schematically in Figure 5; and the actual plot of 
percent resistance versus percent rotation resulting is indicated by 
circles in the figure. It is noted that 61.4% of the total resistance 
will now occur at about 75% of total potentiometer rotation. If this point 
is denoted for the simulator as having a control rod position reading of 
zero, then a discrepancy will exist between the control rod-reactivity 
relationships of the actual reactor and that of the simulator. As seen 
from Figure 5, the simulator rods increase reactivity linearly over their 
entire travel while actual rod insertion increases reactivity linearly 


only over the last 10 cm of travel. This situation must be accepted when 
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the actual reactivity curve cannot be duplicated over its entire range, 
The rate of potentiometer rotation is easily obtainable by timing the . 
actual control rod insertion rate, Both fine and coarse control rods have 


been set to move from the fully retracted to the fully inserted positions 2 
(O to 25 cm.) in about one minute at the fast speed setting, Therefore, / 
the last 10 cm of actual rod travel will be completed іп 0,40 min, From 
Figure 5 the last 10,cm of actual rod travel commences at 85% of total 
potentiometer rotation or, alternately, over the last 0,15 revolution, 
Hence, the reactivity potentiometer drive rate should be 0,15 revs/0,40 
nmin = „37 rpm, at the proper rate, Over this straight line portion, 


reactivity change at fast speed is found to be 
2 4 
837 к? /em x - cm/sec = 349 у? /зес 


To solve the second problem, & means of establishing a variable zero 
or critical reference level must be provided for the reactivity potentio- 
meter, The reactor functions between certain definite reactivity limits 
for any specified core temperature, In other words, some scheme is re- 
quired to duplicate the normal operating ranges to be expected from 
maximum suberitical conditions to maximum supercritical conditions, For 
example, at 23°C, the maximum positive reactivity attained with all rods 
fully inserted is only about „0020, At the same time, with all rods fully 
retracted, the reactor has a negative reactivity of about „035. This 
one-sided variation is representative of the type of reactivity behavior 
| encountered under typical operating conditions, It should be noted here 
that the above values were determined with the so-called "glory hole" empty, 


The glory-hole (Figure 2) is a small access tube about one inch in diameter 


27 





which passes through the core center, and into which additional fuel, 
moderator, reflector material, samples and the like can be inserted to 
yield changes in reactivity. As an example, by insertion of the proper 
amounts of moderator, reflector, and lead shielding in the glory-hole, 

the core is made essentially homogeneous; and an additional reactivity 
contribution of about 2600 ре can be realized. Moreover, an Tn 
linear temperature coefficient of reactivity, based on water shield tem- 
perature, of about -2.75 x 107^ рег 9С exists for the AGN-201. This rather 
significant effect on reactivity must also be taken into consideration. 
Figure 6 represents the scheme devised to yield the proper operating range 
of reactivity from a relatively large negative value to a small positive 
value and, in addition, to permit control over maximum positive reactivity 
available. A permanently imposed offbalance in the values of R, and Ry 
will electrically offset the mechanical mid-position as indicated. The 
values given in the figure are those based on the reactivity conditions 

at 23°C and employing a homogeneous core, The critical or zero position 

is now located at point C. In addition, it is noted that C can be further 
controlled by movement of the maximum reactivity control and temperature 
coefficient (MRCTC) potentiometer. If, for example, the temperature in- 
creases (or alternately, if a sample with a large thermal neutron absorption 
cross-section is placed in the glory hole), the MRCTC potentiometer is moved 
as indicated in the figure. This increases the resistance in loop 2 and 
decreases the resistance in loop 1, causing an unbalance in the feedback 
circuits, thereby producing a decrease in reactivity corresponding to the 
amount of resistance change applied. If, therefore, it is desired to hold 
the reactivity at zero, the reactivity potentiometer must be moved in such 


a direction as to balance the resistances in the two loops. In this case, 
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it would require movement of the arms as indicated. It can readily be seen 

that if the temperature coefficient is offset far enough, the maximum 

positive reactivity available can be reduced until the reactor may no longer 

be brought critical even though all rods are inserted. This corresponds 

to the point marked B in the figure. The original electrical center position 

C, therefore, effectively walks up and down the reactivity potentiometer 

depending on the setting of the MRCTC potentiometer. Of course, maximum 

positive reactivity may likewise be increased to prompt critical conditions 

and above if desired. One of the solutions of the transient eauations for a | 
step change in reactivity is the so-called Inhour equation. The plot of this ነ 
equation using the value of neutron lifetime applicable to the AGN-201 is m у” 


given in Figure 7. Calibration of the MRCTC potentiometer can be made by a | 
ы” 
\ 


when all rods are fully inserted. The corresponding reactivity may then be Б 
ت‎ ЕЕ сия —— а መላ Т... -- መም سوت‎ 3 — —. \ 
picked off an appropriate plot of the Inhour equation using the neutron life- | б 


measuring the stable period (T) for any particular potentiometer setting 


— 


time applicable to the reactor being simulated. Such a calibration curve' is 
given in Figure 8 for the particular MRCTC potentiometer employed here. 
Temperatures shown in the figure are based on experimental data employing a 
homogeneous core and using the negative temperature coefficient of reactivity 
from Table 1. In summary, therefore, the MRCTC potentiometer functions as 
an electrical stop on the upper limit of the reactivity potentiometer which 
provides a scheme to preset any desired maximum reactivity (or alternately, 
any desired core temperature). The effect on the electrical circuit can be 
taken into account quite simply as is indicated in Appendix II. 

The third problem can be solved by y employing stablized operational 
amplifiers on Ap zero drift to give reasonably accurate and repro- 


p solutions at the extremely low starting or source voltages encountered. 
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The latter are of the order of magnitude of a few millivolts if scaling is 
such as to have about four decades of multiplication from source level to 
full power. The maximum output of the amplifiers is lOO volts. A variable 
biasing arrangement is also desirable to ensure zero offset at zero voltage 
input to the amplifier gird. Any one of several suitable methods recom- 
mended by the manufacturer of the amplifiers will, satisfy this requirement, 
Causing additional problems at low operating voltages are the condensers in 
the circuit. It is found that an erroneous flux-time curve results in the 
low voltage region upon the decay of the neutron flux when the reactor 
simulator is scrammed. Analysis of oil, paper, mica, electrolytic, and 
other types of commonly used capacitors reveals a marked dielectric ab- 
sorption effect by which these types hold their charge for a relatively 
long period of time even after direct shorting across the terminals. This 
problem is virtually eliminated by the use of polystyrene or mylar condensers 
—— exhibit characteristically low dielectric absorption properties. 

| The fourth and final problem is significant only in that it is con- 
sidered desirable to simulate reactor responses over the entire operating 
range. This demands computer solutions to the controlling equations at : 
extremely low values of neutron flux which corresponds to the condition in 
the reactor when all rods are retracted, and the existing neutron density is 
dependent on the extraneous neutron source strength. This problem may be 


approached as follows: From equation (15), namely, 


dC ^ 
* dn _ + 2 1 


ас 
it is apparent that under steady-state conditions + and — terms 


drop out leaving 
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ог rearranging, 


У; 


So = 15 (23) 
And from equation (8) and (9), this can be written as 
፲(1 = K ap) 
у= ብር — (24) 


From manufacturer's data or from actual experimental flux determinations, a 
value of n can be calculated at a level where reactivity is known, The 
latter is obtainable from experimentally determined control rod calibration 
data, Employing (23), this will establish a value of 9617 at one particular 
control rod position, From the known value of reactivity at this position, 
K er Тау be determined from (P), Finally, from (24), Spl” is obtained, 


If the value of estimated neutron lifetime as given by the manufacturer is 


then employed, Sj is readily calculated. It is noted from (22) that 


а Ыр Пе ӨҚЕ. ры 
R90 (max) lo 801 + ШЕ с 


or rearranging, 


з E AE ы) 
R(1 + / пах 


Therefore, to obtain the neutron source voltage corresponding to a particular 


source generation rate, the product (So (max) 15 ) need only be determined 


3k 





as accomplished above, and actually the two terms need not be separated 
to determine the voltage, The magnitude of neutron lifetime itself is 
important in the determination of a proper value of the integrating 


capacitor cy. From (22). 


% 


AP 


or, essentially, с. is proportional to neutron lifetime, However, if 

E. order of magnitude of neutron lifetime is known, this is sufficiently 
accurate to yleld excellent results as long as maximum reactivity values 

are limited to ,005 and under, This figure of ,005 is arrived at by evalua- 
ting the stable period corresponding to various values of positive re- 
activity when different neutron lifetimes are considered, The Inhour 
Equation may be employed to make this comparison analytically; and it will 
be found that, at low reactivities, period is independent of neutron life- 
time for a specified reactivity, Furthermore, in view of the fact that the 


4 


stable period corresponding to a neutron lifetime of 10 * sec and a re- 
activity of ,005 is about 2,5 seconds, this condition on the actual 
reactor cannot be duplicated because a period scram is set to function at 
about 3,0 seconds, As a consequence, comparison runs at higher values of 
reactivity to better evaluate neutron lifetime are not feasible, 

Having solved the above four problems, the remaining circuit para- 


meters may be calculated as indicated in Appendix II and the final AGN-201 


Ze 
— A EIER 





simulator can be constructed, Figure 9 is a schematic of the final analog, 


tm. 





and listed in Table 2 below are the various values of circuit components 


referred to in the figure, Also listed are the other factors which appear 





in the scaling factor equation (22). 
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TABLE 2a, CIRCUIT VALUES 


В = О (actual value of pot,) 
В, = 119 Ko 

К. = 129 Ке 

c. = .0007 ህድ R = 102M 
с, = .0207 pf В, с 30M 
وه‎ = .208 yf Ry = 10,6M 
е, = шош Ry = 12,0M 
в: + 2.6 yf К, = 15.4M 
EU «785 pf Re = 102M 
о» = .0039 pf 

В, = № = 1M 

Ry = 10M 

V = 51ዉኛ 


TABLE 2b, ADDITIONAL SCALING FACTOR VALUES 


Port ” 


В = 119K (use for В in equation 22) 


23.0К< (use for P in equation 22) 


/* шах= ,034775@ 23°C (core homogeneous) 


н = .00025 

Ва = .00085 

By = .00241 

в, = «00213 

85 З «00166 

% 5 00025 

ІР = ож 

So = So (pax) = 579 neuts/cm3/sec 


7 





1” ዶሪ 1.0 x 1074 sec FY % near critical 


Figure 10 is a photograph of the completed simulator with linear 
recorder unit attached. Important features and actual circuitry are 
shown in Figure 11 which is a rear view of the simulator console, At 
the left center of Figure 11 is the reactivity potentiometer and drive 
motor assembly; at the upper right portion of the figure is the MRCTC 
potentiometer; and at the bottom portion of the photograph is the entire 
simulator circuit which is constructed to resemble the schematic of 
Figure 9, Along the top of the console are located the operational 
amplifiers and stabiliting amplifiers, 

6, Degree of Simulation Attained 

Figures 12 through 14 represent comparison flux-time plots taken 
from the reactor and the simulator, utilizing the same linear recorder 
so as to ensure an identical time base, The sequence of events describing 
each of the comparison runs is briefly discussed in the succeeding 
paragraphs. 

The sequence of events describing Figure 12 is given in Table 3. 


TABLE 3, SEQUENCE OF EVENTS FOR FIGURE 12 


TIME IN REACTOR SEQUENCE CORRESPONDING 
MINUTES SIMULATOR SEQUENCE 
— Steady flux level; rods & 15 Steady flux level; 
em; glory hole empty; temp, rods @ 10.1 cms; МЕСТО 
21°C, set for proper max, 


reactivity conditions, 


0 Inserted rods @ fast speed. емді rods @ fast 
speed, 
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FIG. ll. REAR VIEW OF SIMULATOR CONSOLE 





0.4 


Pl 


4.25 


7.50 


8,65 


10.15 


12.05 


13,5 


Rods пр; flux permitted to 
increase on 17.5 sec, period 


Flux level „ 76x 10° 

(0.76 x 10 ” ppa scale read- 
ing); rods retracted to 
20.45 cm; flux permitted to 
decay to 1,0 x 10”, 


Adjusted rods for steady 
flux level, 


Quickly inserted moderator, 
graphite, and lead discs іп- 
to glory hole to provide a 
step change of about (+) 

1300 we б of total required 
to make core essentially 
homogeneous, ) 


Retracted rods to 19.93 
cm; allowed Гуж to de- 
cay to 4,0 х 107, 


Inserted remainder of 
poly, graphite, and 
lead in glory hole and 
adjusted rods to steady 
flux level of 4,5 x 10 
г rated power 100 mw, 


Quickly withdrew all 

poly from glory hole to 
provide a step change of 
about .(-) 2600 „o ; аї-- 
lowéd flux level to decay 
to about 0,63 x 10°, 


Inserted rods all the 
way; permitted flux to 
increase to about 4,5 
x 10° and leveled at 
about 4,0 x 106, 


Rods up; flux per- 
mitted to increase 
on 17,5 sec, period, 


Flux level @ 7.6 x 


10° (76 volts scale 
reading) rods re- 

tracted to 22,5 cm; 
flux permitted to, 
decay to 1,0 x 10. 


Adjusted rods for 
Steady flux level, 


Quickly rotated MRCTC 
potentiometer to pro- 
vide corresponding 

(+) 1300 pp step change 
in reactivity, 


Retracted rods to 21,9 cm; 
allowed e. decay to 
about 4,0 x 107, 


Rotated MRCTC potentiometer 
an additional amount corr- 
esponding to about 1300 pp 
to simulate homogeneous core 
conditions and brought to 
Buy flux level of 4,5 x 


шш. 


Quickly reset MRCTC to or- 
iginal setting correspond- 
ing to glory hole empty 
conditions thus providing 
a step change in reactiv- 
ity of about (-) 2600 мо } 


allowed flux to, decay to 
about 0,54 x 10 . 


Inserted rods all the way; 
permitted flux to, increase 
to about 4,5 x lO and ረ 
leveled at about 4.0 х 10, 


Ав indicated in Table 3, the position to which simulator rods are retracted 


in each case is different from that given by the rod position indicator on 
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the reactor console, 


This implies that the reactivity potentiometer con- 


tributes more reactivity per centimeter of movement than do the actual rods, 


By a trial and error procedure of slightly altering the value of the 


paralleling resistance of the reactivity potentiometer, the actual and 


simulator rod positions might be more closely brought into agreement, 


TABLE 4. 


TIME IN 


The sequence of events describing Figure 13 is given in Table 4, 


SEQUENCE OF EVENTS FOR FIGURE 13 


REACTOR SEQUENCE 


CORRESPONDING SIM- 


MINUTES ULATOR SEQUENCE 


0.4 


2.1 


3.0 


Approximately steady 
flux level; rods @ 15 са; 
glory hole empty; temp. 
21°C е 


Inserted rods at fast 
speed, 


Rods up; flux permitted 
to increase on 17.5 
sec, period. 


Retracted rods and 
brought to steady flux 


level of 7.0 x 106, 


Scrammed reactor and 
permitted flux to decay 
to observe effects of 
delayed neutrons, 


Steady flux level; 

rods 8 10.1 cm; MRCTC 
potentiometer set for 
proper maximum reactivity 
conditions. 


Inserted rods at fast 
speed. 


Rods up; flux permitted 
to increase on 17.5 
sec, period. 


Retracted rods and brought 
to steady flux level of 


7.0 x 106, 


Scrammed simulator and 
permitted flux to decay 
to observe effects of 
delayed neutron networks, 


Simulation on scramming the reactor is difficult to reproduce because of 


dashpot action on the rods as they are ejected from the core. This action 


was approximated by manually rotating the reactivity potentiometer to 


the all-rods-out position, 


As a result, the simulator decay curve on 


scramming appears to be steeper than the reactor curve. However, the 
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general nature of the two decay schemes is cuite similar with the flux 


levels of reactor and simulator being 3.4 x 102 and 2.3 x 10? respectively 


four minutes after scramming. 


The sequence of events describing Figure 14 is given in Table 5, 


TABLE 5, SEQUENCE OF EVENTS FOR FIGURE 14 


to provide step change from 


a supercritical condition 
of about (+) 1300 po to 
a subcritical condition of 
about (-) 1300 Bo e 


TIME IN REACTOR SEQUENCE CORRESPONDING SIMULATOR 
MINUTES .:0 .: SEQUENCE 

— Approximately steady Steady flux level rods 
flux levels; rods@ 20.0 @ 22.0 cm; MRCTC potentio- 
cm; glory hole empty; meter set for proper max, 
temp. 21%, react. conditions. 

0 Quickly inserted all Quickly rotated MRCTC potentio- 
poly, graphite, and lead meter to simulate step change 
discs into glory hole to of about 2600 вр bringing 
provide a step change from simulator from a subcritical 
about (-) 1300 pe to condition of about (-) 1300 e 
(+) 1300 pe conditions. to a supercritical condition 
Permitted flux to jncrease of about (+) 1300 p2 . 
to about 0,73 x 10”, Permitted flux level to in- 

| crease to about 0.73 x 106, 
2.5 Quickly removed all poly Quickly rotated MRCTC potentio- 


meter back to original setting 
corresponding to a step change 
in reactivity from about (+) 
1300 цу? to about (-) 


1300 Зай о 


The current output in рр amps. from an ion chamber feeds the linear 
recorder on the reactor; the appropriate scales are labeled in the three 
figures, Corresponding to these are the voltage scales on the simulator 


which are also labeled in the figures, For example, the l x 1078 ኮኮቤ 
scale on the reactor corresponds to the 100 volt scale on the simulator, 


the 1 x 1077 рра corresponds to the 10 volt scale, the 1 x 10710 ppa 
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to the 1 volt scale, and the 1 x I ppa to the 0,1 volt scale. The 


reading.in ppa as well as the reading іп volts, of course, are proportional 


to a certain flux level. The flux level at the center of the core at a 
power level of 100 mw is given as 4.5 x 106. This corresponds to a 
0.45 x 1078 ppa output from the ion chamber or a 45 volt output from 


the simulator; all other values are scaled accordingly. 


7. Conclusions 

As evidenced by the comparison curves of Figures 12 through 14, reason- 
ably accurate simulation of the time dependent behavior of a specific low 
power thermal reactor may be readily accomplished despite the many simplifying 
assumptions made in the development and electrical duplication of the 
governing transient equations. The final simulator with its relatively 
simple control system provides a useful tool with which to observe typical 
reactor responses resulting from the many sequences of operations to which 


a reactor of this type might be subject. 
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APPENDIX 


CIRCUIT ANALYSIS DETERMINATION OF CURRENT La 










P e AO FOR ANY A зы 
—» lg 
= E 

Assume current directions as indicated, Therefore, 

"۱ ۲ з 
Now, let A = proportional distance of P from mid or zero position, 
Then 

ай -E - 0 

1, 7 ፔጉዕኘ2 (1-4) 

uw 0-Е 

a R+P/2 (1.4 


Therefore, 


-E E 


Is а R + p/2 (1 + А) : В + ۳۸/2 41 - А 
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or 


or 


В Ee x 1 
la E E (1 - A) ۲ R + P/2 የ 


۲ 2 A P/2 
f R + P/2 (1-471 ٩ 8 + 2 MR M 


АР 
з ВР + + - = 
R2 i p a A) + (1 a) | + pa га 


de ^ E 
2 
E +R + а-а) 
ለጄ 
m DUE: 
R 0۰ ۰ 02) | 
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APPENDIX II 
SCALING FACTOR CALCULATIONS 


Referring to scaling factor equation, namely 


ee ® 


let Ro = 10 Megohms, This is an arbitrary value and is purposely chosen 





large to permit reasonable values of voltage V which is dependent on the 
very short neutron lifetime of this reactor, 

Next, assume a value of neutron lifetime 1* of 107% вес, Since the reactor 
will not be far from critical at any time, 1% and 15 will be essentially 
equal to each other, Hence, as a first approximation is will also be 
assigned a value of lost sec, 

Now from control rod calibration data, ax can be determined. In view 
of the fact that normal operating temperatures are about 2990, ae 

will be determined at this temperature assuming a homogeneous core (i.e., 
moderator, reflector, and shielding material in the glory hole). The 


operating range of reactivity here is 
(-) .034775 to (+) .004550 


ዘም 0095 መረር can be taken as .034775$ and an imposed unbalance in the 
resistance values of R will prevent any positive reactivity value greater 


than .004550, In other words, when A is (-) ብ” = Penk) 


A relationship must now be established between voltage E and neutron 


population density n, From known data, 
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$ at 100 mw = 4,5 x 10° neuts/cm*/sec 


1.28 x 104 TOK cm/sec 


and у с 
ог у= 2,22 х 10? cm/sec at 23°C, | 
6 
б E х 10 z 3 
Therefore, n = ? = EXE 20.25 neuts/cm 


at rated power level, 100 mv, 

Because the maximum output of the amplifier is 100 volts and because 1% 
might be desirable to directly correlate flux with observed voltage, 45 
volts is made to correspond to 4.5 x 106 neutrons per cm? per sec, The 


constant of proportionality Í is then given by 


J=- E - — = 22 


Next, & sultable value ef potentiometer resistance P is selected so that 
after modifications by an external resistance to & center tap connection, 
a reasonable value of total resistance is still avallable, A 20К рої 
meets this requirement; and after een as indicated in Figure 5, 
an overall resistance of 13,33K results (13.0K measured), It is to be 
noted that this value is not that which is to be substituted for P in 


equation (22), A Popp Must be determined which will be compatible with 


f 
the derivation of the reactivity potentiometer relationships given in 
section 4, As far as the circuit is concerned, for zero or critical set- 
ting of the potentiometer the resistances on each side must be balanced, 


Therefore, since a (-) ,034775 reactivity is required on one side, or in 


terms of resistance 


( 03477 


e - б ነ S 
334775 + 0255) 13 = (.085)*(13) = 11,5 KS2 
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This implies a value of 


Pore Se a TI AREL = 29.0 አፈር 


which should then be utilized in (2?) in place of P, Now, because the 
actual resistance of P is only 13.0 K, then 10K additional resistance is 
required in loop 2 as indicated in Figure 6, or R, = R + 10К.0., Values 
of В, аге readily available and are given in the summary of scaling factor 
values at the end of section 5, Before calculations begin, a value of 

R must also be selected, First, however, since it is on this resistance 
that the individual values in the six Re networks will depend, investiga- 
tion to ensure reasonable magnitudes for the Re delay networks is 
necessary, Based upon this investigation, a final value of R = 119K yields 
suitable magnitudes for all six delay components, Therefore, 

EJ = 119K and R, = 119K + 10K - 129 Ко. , It is also noted here that 22 
5P which yields better than 1¢ linearity between jo and A, Calculations 


can now proceed as follows; 





R(1 + RP)O = 119 x10? (1 є 342) 0.034775 
= EE 
Now, from (22) 
i са й 
чи R(1 + / тах 


29 





+ 2 
10 Б 


J 
D RQ + Р) max 2,25 Ж 10% 


Cr > 579592 а 
Likewise, 
á _ 
R (1 + КР) ax By 84 
Es: R (1 + R/P)/O ax = "2 o: 10^ 
Bs P4 
1 = R E 
By A, sec , megs 9, p 
n. ¿00025 14 102 0.0007 
2 .00085 1,61 30 0,0207 
3 00241 0,455 10,6 0,208 
4 00213 0,154 12,0 0.5447, 
5 00166 0,0315 1524 2,06 
6 00025 0,0125 102 0.785 





An 


= „00755 


Next, the method for deciding upon a suitable value of source voltage, 


V, as discussed in section 5, is as follows: 


۳ n(1-K, np) ы PS 
D O Ж 
1 lo 
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Now, 8:30 ie Power ratioxn 6) Com 
30 x 1079 
ር ን ሥሥ 5... 
100 х 10 
or n 6.08 х 1072 neuts ПАУ ѕес 
@30 рм 
n 
۲ == 6 safeties up position = (-) „009625 
eff 


from control rod calibration data, From this then, 


Keer = .99047 
ог (1-K np) - ‚00953 


Therefore, 


352* = 59 x 10” 


and 515 .5°5 x 107" 


Therefor m 
e ore, $ (S5 max "my 
R 2 - የ 
E 2.22)* (10 x 10€)*(,585 х 1079 
or У = 2 
(2.55 х 10°) 
й = 510 m 


It is to be noted that, in the case of the simulation of the AGN-201, 


Sg and ስ are the same since it is not desired to vary the source 
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strength (voltage); and therefore only the one magnitude determined above 
is needed to represent the effectively constant Ra-Be neutron source in 
the reactor, 

Finally, if 1* is taken as 107% sec, then the value of Sg is readily 


found to be 


-4 
So = === = MD з neuts/cm^/sec 
10 


5€ 
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